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Abstract. The cubical barycentric subdivision sd c (fC) of a cubical complex 
K is introduced as an analogue of the barycentric subdivision of a simplicial 
(-H , complex. Explicit formulas for the short and long cubical h- vector of sd c (K) 

r~{ • are given, in terms of those of K. It is deduced that symmetry and non- 

negativity of these /i-vectors, as well as real rootedness of the short cubical 
h-polynomial, are preserved under cubical barycentric subdivision. The as- 
ymptotic behavior of the short and long cubical h-vectors of successive cubical 
barycentric subdivisions of K is also determined. 



o 
u 

1. Introduction 

The present work is partly motivated by [2]. In that article, Brenti and Welker 
study the transformation of the h- vector of a simplicial complex A under barycentric 
£SJ ■ subdivision. They express the entries of the h- vector of the barycentric subdivision 

of A as nonnegative integer linear combinations of those of the h- vector of A. In 
particular, they show that symmetry and nonnegativity of the /i-vector are pre- 
served under barycentric subdivision. Moreover, they prove that if the h- vector of 
A is nonnegative, then the /i-polynomial of the barycentric subdivision of A has 
only real roots. 

It is natural to inquire whether similar results hold for non-simplicial complexes. 
We study the transformation of the short and the long cubical h- vector of a cubical 
complex K under a natural cubical analogue of simplicial barycentric subdivision, 
which we call cubical barycentric subdivision. The cubical barycentric subdivision 
sd c (K ) is a cubical complex which subdivides K, so that the poset of nonempty faces 
of sd c (K) is isomorphic to the set of closed intervals in the poset of nonempty faces of 
K , partially ordered by inclusion (see Section 2.3 for a precise definition). Our main 
result expresses the entries of the short and long cubical h- vector of sd c (K) explicitly 
as nonnegative linear combinations of the entries of the corresponding /i-vector of 
K (Theorems 13.21 and 14. ip . From these expressions we deduce that symmetry and 
nonnegativity of the short and long cubical h- vector, as well as real rootedness of 
the short cubical /i-polynomial, are preserved under cubical barycentric subdivision. 
We also study the asymptotic behavior of the short and long cubical /i-polynomials 
under successive cubical barycentric subdivisions ( Corollaries 13.81 and 14. 5[) . 
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2. Preliminaries 



This section reviews basic definitions concerning cubical complexes and their /- 
vectors and h- vectors and introduces the concept of cubical barycentric subdivision 
for such complexes. For backround and any undefined terminology on partially 
ordered sets and on polytopes and polyhedral complexes we refer the reader to [4j 
Chapter 3] and to [H [5] , respectively. 

2.1. Cubical complexes. We denote by Cd the standard d-dimensional cube 
[0, if C if Any poly tope which is combinatorially isomorphic to Cd is said to be 
a combinatorial d-cube. A cubical complex is a finite collection K of combinatorial 
cubes in R™ , such that (i) every face of an element of K also belongs to K and (ii) 
the intersection of any two elements of if is a face of both. The elements of K are 
called faces. 

We denote by F(K) the face poset of K, meaning the set of faces of K, partially 
ordered by inclusion. This poset is a meet-semilattice. The empty face is the 
minimum element of F{K) and the vertices of K are its atoms. The maximal 
elements of T(K) are called facets. The dimension of K, denoted by dim (if), is 
defined as the maximum dimension of a face. 

2.2. Face enumeration. Let K be a (d — l)-dimensional cubical complex. We 
denote by fi(K) the number of i-dimcnsional faces of K. The f -vector of K is 
defined as 

f(K) = (f (K)J 1 (K),...Jd-i(K)). 

The polynomial 

is called the f -polynomial of K. The short cubical h-polynomial of K is defined in 
P] by the equation 

d-1 d-1 

(1) hk e \x) = 5> (SC) W^ = ^f j (.K)(2xy(l-x) d - 1 ~ i . 

4=0 j=0 

The vector 

h^\K) = (h^\K)Mr\K) 1 ... 1 h d S l\{K)) 
of coefficients of this polynomial is called the short cubical h-vector of K. The 
polynomials /# (x) and h^ (x) are related by the equations 

(2) h^\x) = (1 _*)*-!/* (j?-- 

and 

(3) 2 d - 1 f K (x) = (x + 2) d - 1 hP 



f K (x) = 5>(if); 



x + 2 

As a result, the entries of f(K) can be expressed in terms of those of h} sc '{K) and 
vice versa by the equations 

'd-l — iV(sc), 



(4) m) = 2- j EL_ .r'w 

i=0 ^ ■'' 
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and 

(5) 



3=0 ^ % ' 



The (long) cubical h-vector 

M c )(a) = (4 c) (a),/4 c) (a),...,4 c) (a)) 

of K is denned in [I] by the recursive formula 

(6) h^ (A) = h\ e) (A) + h<-% (JRQ, for < i < d - 1 
and the initial condition /ig (A) — 2 d ~ 1 . The polynomial 

i=0 

is called the (long) cubical h-polynomial of A". The short and long cubical /i- 
polynomials of A are related by the equation 

(7) (l + x)h$(x) = 2 d - 1 +xhP{x)+2 d - 1 (-x) d+1 x{K), 
where 



d-l 



x(K) = -i+j2(-iyfi(K) = -i+f K (-i) 

is the reduced Euler characteristic of A. The entries of h^ c \K) can be expressed 
in terms of those of hS sc '{K) by the equation 



(8) 



hf\K) = ^{-l) i+ i- l h { * c \K) + (-iy2 d -\ 1 < i < d. 

3=0 

The cubical /i-vectors of A can also be expressed in terms of the simplicial h- 
vectors of the links of the vertices of A (these links are simplicial complexes) ; see 
1, Theorem 9]. 

2.3. Cubical barycentric subdivision. Let A" be a cubical complex. The cubical 
barycentric subdivision sd c (A) of A is the polyhedral complex defined as follows: 
Vertices of sd c (A) are the barycenters of the nonempty faces of A. Furthermore, 
to each closed interval [F, G] in the poset F(K) \ {0} of nonempty faces of A 
corresponds a face of sd c ( A) . This face is equal to the convex hull of the barycenters 
of all elements of [F, G]\ see Figure 1 for an example. 



i «< * 4 i 

1 .. «. 

■ i|i * H 



Figure 1: A cubical complex and its cubical barycentric subdivision 
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The face poset of sd c (^) is isomorphic to the set of closed intervals in the poset 
J~(K) \ {0} of nonempty faces of K, partially ordered by inclusion, with a minimum 
element adjointed. This fact implies that sd c (K) is also a cubical complex; its 
dimension is equal to that of K . 

3. The short cubical /i-vector 

This section studies the short cubical /i-polynomial of sd c (i^) with respect to 
its nonnegativity, symmetry and real rootedness. The asymptotic behavior of the 
short cubical /i-polynomial of the n th iterated cubical barycentric subdivision of 
K, as n approaches infinity, is also explored. Here and in the following section, K 
denotes a (d — l)-dimensional cubical complex and sd c (K) stands for its cubical 
barycentric subdivision. 

To express the short cubical /i-vector of sd c (K ) in terms of that of K, we first 
need to determine the relationship between the corresponding /-vectors. 

Proposition 3.1. The f -vectors of K and sd c (K) are related as follows: 

d-l , .s 

(9) fi(sd c (K)) = ^PUOfQ, i = 0,...,d-l. 

Proof. We recall that there is a one to one correspondence between the set of i- 
dimensional faces of sd c (K) and the set of closed intervals [x, y) of rank i in the face 
poset of K. To count these intervals, we note that there are fj(K) ways to choose 
a face y of K of given dimension j. Since every such face is a combinatorial j'-cube, 
there are 2*(^) ways to choose a face x of y of codimension i. As a result, there 
are 2 l (^jfj(K) intervals [x,y] of rank i in the face poset of K with dim(y) = j. 
Summing over all j >i, we obtain ([9]). □ 

Theorem 3.2. The short cubical h-vectors of K and sd c (K) are related as follows: 

d-l 

(10) ^ (sc) (sd c (X)) = ^i?(d, ? ,j)^ c) (K), t = 0,...,d-l 1 

3=0 

where the coefficients B(d,i,j) are nonnegative rational numbers, determined by 
the generating function 
d-l 

(11) £>(<*, mV = ^(sx + iyix + sf- 1 -! . 

i=0 

Proof. Equations Q can be rewritten as 

(12) fsd e (K)(x) = f K (l + 2x), 

where fic{x) and f s d c (K)( x ) are the /-polynomials of K and sd c (K), respectively. 
Using ([2]) and ©, we can further rewrite (fT2")) as 

(13) ^h^wix) = (x + 3) d -^ (|±i) , 
so that 



tiSwto = i E^ c) w( 3 -+ 1 ) J (-+ 3 )^ 1 ^- 



3=0 
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The result follows by equating the coefficients of x % in the two hand sides of this 
equation. □ 

Recall that a vector (ao,ai, . . . , a n _i) 6 K™ is said to be nonnegative (respec- 
tively, symmetric) if a% > (respectively, a; — a„_i_i) holds for < i < n — 1. 

Corollary 3.3. //if /ias nonnegative short cubical h-vector, then so does sdc(-K'). 

Proof. Equation (fTTj) implies that the coefficients B(d,i,j) are nonnegative for all 
< i,j < d — 1. Thus the statement follows from (fTu| . D 



Corollary 3.4. If K has symmetric short cubical h-vector, then so does sd c (K). 

Proof. Replacing x by 1/x in (jlip and multiplying by x , we find that 

B(d,i,j) = B(d,d-l-i,d-l-j). 

Assume that h {sc \K) is symmetric, so that fc,^ .(if) = h^ (K) holds for < 
j < d - 1. Then 

d-l 



fttl-^BdcCif)) = X)B(d,d-l-*,j)ft$* c) (liO 

i=o 

= £B(d,d-l-i,d-l-/c)/4-L 

fc=0 

= ^i?(d,z,fc)4 sc) = ^ sc) (sd c (if)) 



fc=0 

and hence /i( sc )(sd c (if)) is also symmetric. □ 

Corollary 3.5. The short cubical h-polynomial of K has only real roots if and only 
if the same holds for sd c (if ). 



Proof. This statement follows easily from (|13[) . The details are left to the reader. 

□ 

Example 3.6. For the boundary complex K of the 3-dimensional cube we have 
h£ c (x) = 8(1 + x + x 2 ). This polynomial has positive coefficients and two non-real 
complex roots. By Corollaries l3.3l and l3.5l so does the short cubical ^-polynomial of 
sd c (if). This situation is in contrast with what holds for barycentric subdivisions 
of simplicial complexes (see [21 Theorem 2]). 

We denote by sd™(if) the n th iterated cubical barycentic subdivision of if, i.e. 
sd° c (K) = K and sd"(if) = sd c (sd™ _1 C*0) for n > 1. The short cubical h- 
polynomial of sd"(if) has the following simple expression, in terms of the short 
cubical /i-polynomial of if. 

Proposition 3.7. The short cubical h-polynomials of K and sd"(if) are related as 
follows: 

U4) ^ W W-^ 2 ) hK v (2»-l)x + 2" + lJ- 
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Proof. This follows from (fT3| by induction onn. □ 

The following statement implies that all complex roots of the short cubical h- 
polynomial of sd"(if) converge to — 1 as n — > oo. 

Corollary 3.8. We have 

coefficientwise, as n —> oo. JVi particular, the short cubical h-polynomial ofsd™(K) 
has positive and unimodal coefficients for all large n. 

Proof. The first statement follows from flU]) and the fact that h^ c \l) = 2 d ~ 1 f d -i(K) 
by straightforward computations. The second statement follows from the first and 
well-known unimodality properties of binomial coefficients. □ 

Remark 3.9. We do not know of an example of cubical complex K for which 
h^ sc \K) is nonnegative and h^ sc ^ (sd c (K)) is not unimodal. 

4. The cubical /j-vector 

This section proves results on the cubical ft,- vector of sd c (K) analogous to those 
on the short cubical h- vector in Section [3] 

Theorem 4.1. The cubical h-vectors of K and sd c (K) are related as follows: 

d 

(15) h[ c) (sd c (K)) = J2C(d,iJ)hf(K), 

j=o 
where the coefficients C{d, i,j) are nonnegative rational numbers, determined by the 
generating function 

ei«uw = T^-d + ««v) + a (x + 3) T 3 ~ % +T v " 

*ri. ~^ I + x 2 aj x + 3 — (Sx + l)y 

(16) + 2d _ 1(1 + x) ((* + if- l + (3* + l) d -V). 

Proof. Since ^ c) (iif) = /i[, c) (sd c (i<)) = 2 d -\ equation (jS} is valid for ? = if we 
set C(d, 0, 0) = 1 and C(d, 0, j) = for 1 < j < d. This agrees with ([15]). since the 
right-hand side reduces to the constant polynomial with value 1 for x — 0. Using 
©, © and (JTUJ), we find that for 1 < i < d 



h^(sd c (K)) = ^(-l)'^- 1 4 sc) (sd c (X)) + (-l) l 2 rf - 1 

fc=0 

= ^(-l) i+k - l J2B(d,k,j)h^\K) + (-l)^- 1 

fc=0 j=0 

i-1 d-1 

= ^(-i) i+fe - 1 ^B(^,i)(/ l f(^)+ftSti(^)) + (-i) j 4 c) w 

fc=0 j=0 

d 

3=0 
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where 

i-1 

(17) C(d,i,0) = ^(-l) i+fc - 1 B(d,fc,0) + (-ir 

fc=0 

and 

i-l 

(18) C(d,i,i) = ^(-l^-^B^fe.jO + B^fe.i-l)) 

fc=0 

for 1 < j < d, under the convention that B(d,k,d) = for all k. Using (fTTj) . we 
compute that for 1 < j < d — 1 

d d i— 1 

^w^v = ^^(-lr+^WA^ + s^j-i))^ 

i=0 1=0 fc=0 



d ~ 1 1 _ (_ v \d-k 

j 1 + i 

d-\ d-\ 



fe=0 



1- (x J2 B(d, k, j)x k + (-x) d+1 J2 B(d, k, j)(-l) k 

V fc=0 fe=0 

d-1 d-1 \ 

xJ2 B (d,k,J - l)^ fe + N W E B (^.J - 1)(-1)* 

fe=0 fe=0 / 

(s(3» + l) J (x + S)^ 1 ^' + (-.x) d+1 (-2) J 2 rf - 1 - J 



2 d -!(l+x) 

+ x(3x + l) J '- 1 (a; + 3) d - J + (-.x) d+1 (-2)J- 1 2 d -J) 



and hence 

d 

(19) 2C(d,t,jV = 



a;(3a; + iy- 1 (a; + 3) d - 1 -^ 



2d-3 

i=0 



Similar computations yield 



ii 



1 /g(38 + l)^ ^ , +1 



(20) E^^^Tb C^x" +1 ) 

and 

d 

(21) £<W,d)*' = M 2d _ 3 

i=0 V 

From these equations we can infer that the C(d, i, j) are nonnegative rational num- 
bers for all < i,j < d. Multiplying the equations (f!9|). ([20]) and (J2T|) by y J for 
l<i<rf— 1, i = and j = d, respectively, and summing over all j results in (IT51) . 

D 

Corollary 4.2. If K has nonnegative cubical h-vector, then so does sd c (K). □ 
Corollary 4.3. If K has symmetric cubical h-vector, then so does sd c (K). 



8 CHRISTINA SAVVIDOU 

Proof. As in the proof of Corollary 13.41 it suffices to show that C(d, d ~ i,d — j) = 
C(d,i,j) for all < i,j < d. This follows from (|16[) by replacing x and y by l/x 
and 1/y, respectively, and multiplying by x d y d . □ 

Remark 4.4. We do not know of an example of cubical complex K with nonneg- 
ative cubical h- vector for which %V (K ^(x) is not real-rooted. 

Corollary 4.5. For d > 2, we have 

(22) aiiy^W*) "> /*-!(*)* (S + 1 ) 41 " 2 

coefficientwise, as n — > oo. In particular, if (-l)x(K) > 0, then the cubical 
h-polynomial ofsd"(K) has nonnegative and unimodal coefficients for all large n. 

Proof. Since sd c (K ) is a subdivision of K , we have x(sd c (if)) = x{K) (this equality 
also follows from (TH21) by setting a; = —1). Thus, by applying (J?) to sd™(if), we get 

(23) (l + 4 (Jf) W = 2 d - l +xh^} (K) {x)+2 d -\-x) d ^ X {K). 

Note that if (-l) d - 1 x( J ftT) > 0, then ^^(sd"^)) = {~2) d - l x{K) > 0. The result 
follows from {23]) and Corollary [3JH1 D 
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